Method and apparatus for reducing crosstalk in a DWDM transmission system

ABSTRACT

A method and apparatus is provided for transmitting a WDM optical signal. The method begins by modulating a plurality of optical channels that are each located at a different wavelength from one another with a respective one of a plurality of information-bearing electrical signals The channel spacings between adjacent ones of the plurality of optical channels are selected such that the optical channels give rise to selected FWM terms that lie outside a channel bandwidth of the optical channels. Each of the modulated optical channels is multiplexed to form a WDM optical signal. The WDM optical signal is then forwarded onto an optical transmission path.

FIELD OF THE INVENTION

The invention relates generally to the transmission of wavelength division multiplexed (WDM) optical signals, and more particularly to a method and apparatus for reducing nonlinear signal distortions induced by four wave mixing (FWM) interactions between the optical channels.

BACKGROUND OF THE INVENTION

MSOs (multi-service operators) of cable television (CATV) systems provide several services to end users through a fiber optic network, with the final connection to the user often being through a coaxial connection. The services provided by the MSO typically include broadcast analog video and narrow cast digital services, such as data, VoIP, subscription, pay per view and video on demand (VOD) services. The services are generally allocated to a portion of an optical channel, which typically has approximately 1 GHz of available bandwidth. While the bandwidth of a channel is generally constrained by the optical network (the optical network transmitters, optical fiber, channel filter bandwidths and coaxial connection), the number of users connected to the optical network continues to increase, which requires increased demand for bandwidth for the desired services.

In recent years wavelength division multiplexed (WDM) optical transmission systems have been increasingly deployed in optical networks. These include coarse wavelength division multiplexed (CWDM) and dense wavelength division multiplexed (DWDM) systems. Whether a system is considered to be CWDM or DWDM simply depends upon the optical frequency spacing of the channels utilized in the system. Although WDM optical transmission systems have increased the speed and capacity of optical networks, the performance of such systems is limited by various factors such as chromatic dispersion and the fiber nonlinearity, which can cause pulse shape changes in the case of baseband digital signals and distortions in case of analog signals. These impairments degrade the quality of the optically transmitted information. Fiber nonlinearities, for example, can give rise to the generation of undesirable optical signals whose wavelengths fall within the filter passbands of the various channels within the system. If the power levels of these undesirable signals are sufficiently large they will corrupt the information being transmitted by the optical signals that are impacted by these impairments.

One common nonlinear impairment in an optical fiber communication system having multiple wavelengths is four wave mixing (FWM). The FWM process is also known by other names such as parametric interactions. Parametric interactions, or FWM, occur when the photons comprising the desired optical signals of a system elastically scatter in a silica fiber (and other materials). The material comprising the optical fiber acts as catalyst for elastic scattering. During the course of elastic scattering events incident photons interact with the propagation medium (the optical fiber material) through collisions with the material. The scattering event is elastic because none of the optical energy transported by the incident photons is lost to the fiber itself. This is unlike Raman scattering, which is an inelastic scattering process with some of the optical energy being lost to the transmission medium. Instead, the incident photons' energy is instantaneously converted back into optical energy. The scattered photons' energy may or may not be at the original frequencies of the incident photons that initiated the elastic scattering event. Both energy and momentum are conserved during the elastic scattering event:

ω₁+

ω₂=

ω₃+

ω₄ and

{right arrow over (k)} ₁ +

{right arrow over (k)} ₂ =

{right arrow over (k)} ₃ +

{right arrow over (k)} ₄, where ω_(i) is the frequency of the i^(th) signal and k_(i) is its associated wavevector,

is Planck's constant divided by 2π. As a result of the scattering event two photons are annihilated and two other photons are created. If the energy relationship is rearranged it becomes clear how three co-propagating waves interact to create a fourth signal: ω₁+ω₂−ω₃=ω₄. The fourth signal at ω₄ is a new signal that was not originally present in the system at the transmitter launch site. The signal at ω₄ is created as a result of the inelastic scattering events (four wave mixing processes) occurring within the optical fiber. The fourth signal at ω₄ is the undesirable nonlinear impairment. If ω₄ is too close in frequency to one of the original signals (i.e. one of the DWDM channels) it can potentially interfere with and corrupt the information being transmitted by the desired signal in its vicinity. The consequence of momentum conservation is that the four wave mixing terms are largest under conditions when the dispersion is smallest amongst the interacting optical signals. This situation arises when the optical DWDM signals are located near the zero dispersion point of the optical fiber (i.e. in the vicinity of 1311 nm for ordinary single mode optical fiber) or when the signals are closely spaced together (i.e. when the signals are spaced at distances of 100 GHz or less). In the case of a near zero dispersion system the optical pump and signal waves are propagating at nearly identical group velocities through the media. The zero dispersion wavelength of a transmission media refers to the wavelength at which an optical signal will have no change in (inverse) group velocity with respect to changes in its optical frequency. The zero dispersion wavelength differs for different transmission media. The dispersion will generally increase as the wavelength difference between the optical signal and the zero dispersion wavelength increases. However if the signals are closely spaced in frequency the difference will not be great enough to disrupt the conservation of momentum requirement to any great degree.

As noted above, four-wave mixing is a nonlinear process in which the interaction of three fields leads to the generation of a fourth field. That is, FWM is the induced combination of three wavelengths to produce one or more new wavelengths. Two of the three combining wavelengths can be degenerate such that two wavelengths are combined to produce one or more new wavelengths. Optical power is removed from the combining wavelengths and transferred to the new wavelengths in the FWM process. FWM is especially problematic in optical communications if the new wavelengths produced by the FWM process overlap the assigned wavelengths of existing DWDM channels because it is difficult to distinguish between the legitimate optical data signals at these existing WDM channels and the error signal superimposed thereon as the result of FWM. Therefore, it is important to suppress FWM.

FIG. 1 illustrates the four-wave mixing process. If three signals at frequencies f1, f2, and f3 propagate through a single mode fiber, light at a frequency f4=f1+f2−f3 will be generated. If any two of the frequencies are identical the FWM process is referred to partially degenerate and if all three frequencies are identical the FWM process is referred to as totally degenerate. More generally, in a system with N channels, three signals at frequencies fi, fj and fk produce signals at frequencies fi±fj±fk, where i, j and k vary from 1 to N. This results in N²(N−1)/2 interfering signals. FIG. 2 shows an example with three channels, which produce 12 interfering terms.

Optical WDM networks typically allocate a portion of the spectrum about a center frequency of the nominal channel wavelength for signal transmission. Two different optical wavebands or windows are commonly utilized, one centered around a wavelength of about 1310 nm and another centered around a wavelength of about 1550 nm. Signals operating in the 1310 nm window typically occupy wavelengths between about 1280 nm and 1330 nm, whereas signals operating in the 1550 nm window typically occupy wavelengths between about 1530 and 1565 nm. In dense wavelength division multiplexing (DWDM) systems, channel spacings of less than 1 nm are typically used with wavelength bands centering usually centered around 1550 nm. In the 1550 nm window, the International Telecommunication Union (ITU) specifies a DWDM standard using 100 GHz or 200 GHz spacing between signal wavelengths, which is equivalent to 0.8 nm and 0.4 nm, respectively.

The widely used optical wavelengths in a typical CATV application, e.g. those in the 1310 nm window, exhibit little relative dispersion between adjacent ITU frequencies, and hence are particularly affected by the FWM effect when used in a DWDM system. For this reason CATV MSOs have been generally prevented from using a DWDM approach to increase bandwidth because of the degradation of the signal quality from FWM interactions. As a result, in order to meet the increased demand for additional bandwidth, CATV MSOs may be required to install more optical fiber to carry additional channels, and then segment their subscriber base between the newly installed optical fiber and the existing fiber. However, this approach requires a significant capital investment for the MSOs and often involves the negotiation of additional access rights to install the optical fiber. Alternatively, the CATV MSOs may use other wavelengths which are less affected by FWM such as those in the 1550 nm window. However use of these wavelengths require more expensive optical components, e.g. lasers and nodes, significant changes to their existing optical network and significant capital investment as well. Hence it remains desirable to operate DWDM systems in the 1310 nm window for a number of reasons. This is primarily because systems operating in the 1310 nm window are generally simpler and less expensive to implement than systems operating in the 1550 nm window. For example, laser transmitters operating in the 1310 nm window are significantly less expensive than lasers operating in the 1550 nm window. Moreover, because of the higher dispersion in the 1510 nm window, additional components such as dispersion compensators are often needed when operating in the 1550 nm window.

Accordingly, it is desirable to reduce the levels of FWM distortions that arise among the individual channels in a DWDM optical system. This is particularly true in the case of a system utilizing optical channels that are located near the zero dispersion wavelength of the transmission medium.

SUMMARY OF THE INVENTION

In accordance with the present invention, a method and apparatus is provided for transmitting a WDM optical signal. The method begins by modulating a plurality of optical channels that are each located at a different wavelength from one another with a respective one of a plurality of information-bearing electrical signals The channel spacings between adjacent ones of the plurality of optical channels are selected such that the optical channels give rise to selected FWM terms that lie outside a channel bandwidth of the optical channels. Each of the modulated optical channels are multiplexed to form a WDM optical signal. The WDM optical signal is then forwarded onto an optical transmission path.

In accordance with one aspect of the invention, the selected FWM terms may include completely non-degenerate FWM products.

In accordance with another aspect of the invention, the channel spacings between each channel i and its adjacent channel may be equal to a constant frequency offset plus a frequency offset δf_(i), wherein the frequency offset δf_(i) is

${{\delta\; f_{i}} \equiv {a_{i}\frac{B}{2}}},$ where B/2 is a minimum acceptable frequency separation and α_(i) is a scaling factor that satisfies the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1

where 1≦i+j−m≦n m≠i≠j, n is the number of optical channels, and i and j are indices denoting any of the optical channels between channel 1 and channel n, which, along with a given one of the optical channels m, produce the selected FWM terms.

In accordance with another aspect of the invention, α₁ may be equal to zero.

In accordance with another aspect of the invention, the scaling factor α_(i) may satisfy only unique ones of the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1.

In accordance with another aspect of the invention, all values of the scaling factor α_(i) may be equal to zero or a positive integer.

In accordance with another aspect of the invention, B may be equal to the channel bandwidth of the optical channels.

In accordance with another aspect of the invention, the optical transmission path may be located in a HFC network.

In accordance with another aspect of the invention, the optical transmission path may be located in a CATV transmission network.

In accordance with another aspect of the invention, the optical transmission path may be located in a PON.

In accordance with another aspect of the invention, the optical channels may be located at wavelengths at or near a zero dispersion wavelength of the transmission path.

In accordance with another aspect of the invention, the optical channels may be centered about a wavelength of 1310 nm.

In accordance with another aspect of the invention, the optical channels may be centered about a wavelength of 1550 nm.

In accordance with another aspect of the invention, a method and apparatus is provided for determining the frequencies of optical channels in a WDM system. The method begins by selecting an initial set of channel spacings between at least three optical channels to be employed in the WDM system. The initial set of channel spacings are adjusted so that the optical channels give rise to completely non-degenerate FWM products that lie outside a channel bandwidth of each of the optical channels.

In accordance with another aspect of the invention, a WDM optical transmitter is provided. The transmitter includes a plurality of optical sources for generating optical channels located at different wavelengths. The transmitter also includes a plurality of optical modulators each having a first input for receiving a respective one of a plurality of information-bearing electrical signals and an output coupled to an input of a respective one of the plurality of optical sources to thereby provide a plurality of modulated optical channels. The channel spacings between adjacent ones of the plurality of optical channels are selected such that the optical channels give rise to selected FWM terms that lie outside a channel bandwidth of the optical channels. A multiplexer is coupled to the plurality of optical sources to receive and combine the modulated optical channels to produce a multiplexed optical signal.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates the four-wave mixing process in which signals at frequencies f1, f2, and f3 generate a FWM term or product at a frequency f4.

FIG. 2 shows an example of the FWM process in which three optical channels produce 12 interfering terms, including degenerate interfering terms.

FIG. 3 shows the system architecture of a passive optical network (PON).

FIG. 4 shows a simplified block diagram of a WDM transmitter that may be employed in a WDM system such as the PON depicted in FIG. 3.

FIGS. 5 a and 5 b are spectral plots showing examples of a set of four optical channels in the 1310 nm window with irregular channel spacings to reduce crosstalk from FWM and Raman crosstalk.

FIG. 6 is a spectral plot showing another example of a set of four optical channels in the 1310 nm window with irregular channel spacings to reduce crosstalk from FWM and Raman crosstalk.

FIG. 7 is a plot showing a set of four optical channels in the 1310 nm window with a regular channel spacing.

FIGS. 8-11 are showing additional examples of sets of four optical channels in the 1310 nm window with irregular channel spacings to reduce crosstalk from FWM and Raman crosstalk.

FIG. 12 is a flowchart showing one particular procedure for selecting the channel spacings in a WDM system having n optical carrier channels so that FWM is reduced.

DETAILED DESCRIPTION

A method and system is described herein for the suppression of nonlinear signal distortions associated with the propagation of a dense wavelength-division multiplexed (DWDM) optical signal. The particular nonlinear signal distortions that are addressed include four-wave mixing (FWM). Although the concepts described herein will be described within the context of a cable television (CATV) network, the techniques may be advantageously employed within various other networks or systems wherein the suppression of nonlinear effects associated with the propagation of an optical signal is desirable. Such networks include, without limitation, various all-optical networks, hybrid fiber-coax (HFC) networks, fiber deep networks, and networks utilizing a passive architecture, which are often referred to as Passive Optical Networks (PONs).

In typical HFC architectures, a broadcast signal is split at optical hubs and then sent to different nodes together with narrowcast signals. At the fiber node the optical signal is converted into an electronic signal, and carried over multiple coax buses for distribution throughout a neighborhood. On the other hand, in a PON architecture, fibers carry signals from an optical line terminator (OLT) to optical splitters, and further to optical-network units (ONUs), where optical-to-electronic conversion takes place. In the case of PON architectures, broadcast and narrowcast signals are sent in a similar manner. Both HFC and PONs generally carry the same downstream signals to multiple customers. In both networks multiple paths are typically used beyond the first node. Use of WDM in a fiber to a node allows further and more node segmentations.

As previously mentioned, crosstalk arising from four-wave mixing is particularly problematic at low dispersions. This result occurs because FWM depends on the phase relationship between the interacting signals (momentum conservation). If all the interfering signals travel with the same group velocity, as would be the case if there were no dispersion, the effect of FWM is reinforced. On the other hand, with dispersion present, the different signals travel with different group velocities. Thus, with increasing dispersion, the different wavelengths alternatively overlap in and out of phase and the net effect is to reduce the mixing efficiency.

Commonly available optical fiber operating in the 1310 nm window has a zero-dispersion wavelength at about 1310 nm. That is, signals traveling in such fibers at wavelengths in the vicinity of 1310 nm experience little or no dispersion, thus increasing the adverse effects of FWM. Nevertheless, DWDM systems operating in the 1310 nm window remain desirable for the previously enumerated reasons. While the concepts and techniques describing herein are thus particularly applicable to networks that operate in the 1310 nm window. However, they are also applicable to networks that operate in other wavebands such as the 1500 nm, window, for example.

FIG. 3 shows the architecture of a PON in its most generalized form. The PON 100 includes a hub 102, remote nodes 104 that are deployed in the field, and network interface units (NIUs) 106. The hub 102, remote nodes 104 and NIUs 106 are in communication with one another over optical fiber links. If the PON 100 is used in a CATV network, hub 102 is a headend, whereas in a telecommunications network, hub 102 is a central office. The NIUs 106 may be terminal equipment located on the customer premises or they may serve multiple customers, in which case the NIUs 106 simply provide another level in the network hierarchy below the remote nodes. A primary advantage of a BPON is its reliability, ease of maintenance and the fact that the field-deployed network does not need to be powered. Accordingly, PONs are often used as access networks by cable TV in CATV networks and telecommunications providers for the purpose of distributing their services from their facility to the customer premises (e.g., a home or business). One example of a PON is sometimes referred to as a Broadband PON (BPON), which is a combination of a PON with wavelength division multiplexing (WDM) for downstream and upstream signals. The WDM techniques can be used for downstream signals and allow different optical wavelengths to support broadcast and multiple narrowcast (dedicated for each wavelength) transmissions on the fiber employed in the BPON.

FIG. 4 shows a simplified block diagram of one example of a WDM transmitter 200 that may be employed in a WDM or DWDM system. For example, the transmitter 200 may be employed in the hub or headend 102 of the network shown in FIG. 3. In operation, data or other information-bearing signals are applied to the inputs of modulators 210 ₁, 210 ₂, 210 ₃, and 210 ₄. The modulators 210 ₁, 210 ₂, 210 ₃, and 210 ₄, in turn, drive lasers 212 ₁, 212 ₂, 212 ₃, and 212 ₄, respectively. Lasers 212 ₁, 212 ₂, 212 ₃, and 212 ₄ generate data modulated optical channels at signal wavelengths λ₁, λ₂, λ₃ and λ₄, respectively, where λ₄>λ₃>λ₂>λ₁. A wavelength division multiplexer (WDM) 214 receives the optical channels and combines them to form a WDM optical signal that is then forwarded onto a single optical transmission path 240. While the WDM transmitter shown in FIG. 4 multiplexes four optical channels onto a single path, those of ordinary skill in the art will recognize that any number of optical channels may be multiplexed in this manner. As shown, in the context of a CATV network, optical channels λ₁, λ₂, λ₃ and λ₄ may be broadcast signals that all contain the same video signals, plus narrowcast signals that are different for different optical channels of λ₁, λ₂, λ₃ and λ₄. Narrowcast signals that are RF frequency multiplexed into broadcast channels are normally much lower in amplitude than broadcast video signals. The arrangement of sending the same broadcast signal and different narrowcast signals over multiple wavelengths (WDM) is a means of providing more segmentation in CATV networks. If both broadcast and narrowcast signals are being transmitted, then transmitter 200 also includes an RF splitter 216 that splits the broadcast signal among the lasers 212 ₁, 212 ₂, 212 ₃, and 212 ₄. The lasers 212 also each receive a different narrowcast signal.

The deleterious effects of FWM can be reduced by irregularly distributing the signal wavelengths within a selected waveband (e.g., the 1310 nm window). In particular, by using carefully selected irregular or nonuniform channel spacings instead of a uniform channel spacing, the FWM products that are generated fall outside the bandwidths of the individual signals. This approach to the reduction of FWM can be demonstrating using four signals or channels that have a nominal channel separation of 200 GHz. Of course, the techniques described herein are more generally applicable to any number of channels with many different nominal channel separations. After illustrating the process for channels, a generalization of the process will be presented that can be used in connection with any number of channels.

For purposes of this analysis, the four signals can be represented by S(t)=cos [2πf ₀ t]+cos [2π(f ₀ +Δf+δ ₁)t]+cos [2π(f ₀+2Δf+δ ₂)t]+cos [2π(f ₀+4Δf+δ ₃)t] where f₀ is the optical carrier frequency of the longest wavelength signal, Δf is a standard DWDM channel separation (for example, in the current analysis Δf=200 GHz in accordance with ITU DWDM standard). The frequencies δ₁, δ₂, and δ₃ are frequency offsets so that the optical carriers no longer fall on frequencies separated exactly by multiples of Δf, also (δ₁, δ₂, δ₃)<Δf.

The resulting four-wave mixing signal is given by (disregarding the absolute amplitudes as only the frequencies are important for this analysis): S _(FWM)(t)={ cos [2πf ₀ t]+cos [2π(f ₀ +Δf+δ ₁)t]+cos [2π(f ₀+2Δf+δ ₂)t]+cos [2π(f ₀+4Δf+δ ₃)

For a four channel DWDM system operating in the 1310 nm window with the signals defined as above, it turns out there are 48 four-wave mixing terms which can be found by repeatedly applying the appropriate trigonometric identities to the cubed signal S_(FWM)(t). Twenty of these terms fall around the third harmonics of the optical carriers and are therefore not relevant. Of the remaining 28 terms ten are potentially problematic.

Initially channel separations of 200 GHz (f2−f1), 200 GHz (f3−f2), and 400 GHz (f4−f3) were selected with δ₁=δ₂=δ₃=0. This frequency plan does not preclude the generation of frequency terms created by the totally degenerate four-wave mixing process. However, as stated above, there are ten, out of the 28 relevant four-wave mixing terms, which may still cause interference problems with the optical carriers. These ten terms are generated from either the totally degenerate cases or the partially degenerate cases.

Of these ten remaining terms four are generated as a result of the totally degenerate four wave mixing processes. In the case of four optical channels the totally degenerate four-wave mixing processes generates four terms; one falling exactly on each of the four optical carrier frequencies. These terms can never be placed outside of the signal bandwidths. Such terms are always present in any optical transmission situation and are usually not problematic at low to moderate power levels.

In order to ensure that the remaining six terms will fall outside of the DWDM channel filter bandwidths (e.g., ±25 GHz in this example) the optical carrier frequencies must be slightly adjusted from the initially selected values. This is done by determining the appropriate values for the frequency offsets δ₁, δ₂, and δ₃. The six remaining four-wave mixing terms can be placed at least 50 GHz away from the nearest optical carriers by setting δ₁ and δ₂ each equal to 50 GHz (i.e. increasing f2 and f3 each by 50 GHz) This places the terms well outside of the DWDM channel filter bandwidths.

Table 1 lists the relevant frequency expressions involving the optical channels f1, f2, f3 and f4 and δ₁, δ₂, and δ₃. The fourth column of the table shows the displacement between the FWM terms f_(FWM) and the optical carrier frequency fc in terms of the frequency offsets δ₁, δ₂, and δ₃. Two numerical examples are given in the last two columns. In the first example δ₁=50 GHz, δ₂=50 GHz, and δ₃=0 GHz and in the second example δ₁=25 GHz, δ₂=25 GHz, and δ₃=0 GHz. Table 2 lists the optical channel frequency separations for each of these two examples presented in Table 1. Table 3 lists a numerical illustration in the 1310 nm window for each of these two examples. In these examples the carrier signals f1, f2, f3 and f4 are chosen to be 228.745 THz, 228.945 THz, 229.145 THz and 229.545 THz, respectively.

TABLE 1 Relevant frequency expressions including the relationships of the frequency offsets δ₁, δ₂, and δ₃ to the frequency displacement of the four-wave mixing terms relative to each of the optical carriers. FWM Displacement FWM Frequencies Displacement FWM from the Frequencies Displacement Optical from the Closest Four Frequency Carrier Optical Carrier Wave Mixing from (ƒ_(FWM) − ƒ_(C)) (ƒ_(FWM) − ƒ_(C)) Frequencies to the with with Optical Carrier Each Optical Optical δ₁ = 50, δ₂ = 50 δ₁ = 25, δ₂ = 25, Optical Frequencies Carrier Carrier δ₃ = 0 δ₃ = 0 Channel (ƒ_(C)) (ƒ_(FWM)) (ƒ_(FWM) − ƒ_(C)) (GHz) (GHz) f1 ƒ₀ ƒ₀ 0 0 0 ƒ₀ + 2δ₁ − δ₂ 2δ₁ − δ₂ 50 25 ƒ₀ + 2δ₂ − δ₃ 2δ₂ − δ₃ 100 50 f2 ƒ₀ + Δƒ + δ₁ ƒ₀ + Δƒ + δ₁ 0 0 0 ƒ₀ + Δƒ + δ₂ − δ₁ δ₂ − 2δ₁ −50 −25 f3 ƒ₀ + 2Δƒ + δ₂ ƒ₀ + 2Δƒ + δ₂ 0 0 0 ƒ₀ + 2Δƒ + δ₃ − δ₂ δ₃ − 2δ₂ −100 −50 ƒ₀ + 2Δƒ + 2δ₁ 2δ₁ − δ₂ 50 25 f3 ƒ₀ + 4Δƒ + δ₃ ƒ₀ + 4Δƒ + δ₃ 0 0 0 ƒ₀ + 4Δƒ + 2δ₂ 2δ₂ − δ₃ 100 50

TABLE 2 The optical channel frequency separations for the two examples given in table 1. Optical Channel Optical Channel Separations Separations with with Optical Δƒ = 200 and Δƒ = 200 and Optical Carrier δ₁ = 50, δ₂ = 50, δ₃ = 0 δ₁ = 25, δ₂ = 25, Channel Frequencies (GHz) δ₃ = 0 (GHz) f1 ƒ₀ f2 ƒ₀ + Δƒ + δ₁ f2 − f1 = 250 f2 − f1 = 225 f3 ƒ₀ + 2 Δƒ + δ₂ f3 − f2 = 200 f3 − f2 = 200 f4 ƒ₀ + 4 Δƒ + δ₃ f4 − f3 = 350 f4 − f3 = 375

TABLE 3 Original optical frequencies chosen for the 1310 DWM and the updated optical frequencies including the offsets. Original Optical Carrier Frequencies Adjusted Optical Carrier (THz) Frequencies (THz) ƒ₁ = 228.745 (1310.598 nm) ƒ₁ = 228.745 (1310.598 nm) ƒ₂ = ƒ₁ + 0.2 = 228.945 (1309.453 nm) ƒ₂ = ƒ₁ + 0.25 = 228.995, (1309.167 nm): δ₁ = 0.05 ƒ₃ = ƒ₁ + 0.4 = 229.145 (1308.310 nm) ƒ₃ = ƒ₁ + 0.45 = 229.195, (1308.025 nm): δ₂ = 0.05 ƒ₄ = ƒ₁ + 0.8 = 229.545 (1306.030 nm) ƒ₄ = ƒ₁ + 0.8 = 229.545, (1306.030 nm): δ₃ = 0

The examples shown in Tables 2 and 3 are plotted in FIGS. 5 a and 6. FIG. 5 a is a spectral plot showing the frequency placement of the four carrier signals along with the four-wave mixing terms with Δf=200 GHz, δ₁=50 GHz, δ₂=50 GHz, and δ₃=0 GHz. FIG. 5 b, discussed below, is the same as FIG. 5 a except that δ₃=150 GHz. In these two examples the closest FWM terms (other than the totally degenerate terms) are at least 50 GHz away from the carriers. FIG. 6 is an spectral plot showing the frequency placement of the four carrier signals along with the four-wave mixing terms with Δf=200 GHz, δ₁=25 GHz, δ₂=25 GHz, and δ₃=0 GHz. In this case there are three FWM signals that are located 25 GHz away from the optical carriers, which is just at the edge of the DWDM channel filter bandwidths.

One advantage of using DWDM channel spacings such as shown in FIGS. 5 and 6 is that in addition to reducing crosstalk arising from FWM, Raman crosstalk is also substantially reduced. In contrast, in a CWDM system, for example, the Raman gain is much higher. This is because the DWDM channels are located within a bandwidth of no more than about 1 THz.

Additional sets of channels along with their respective FWM terms are presented in the following examples of FIGS. 7-11. In FIG. 7 the channels have a uniform channel spacing Δf=200 GHz and an optical channel filter bandwidth=50 GHz. The shortest channel wavelength is λ_(ref)=1310.000 nm. This plot shows the locations of 10 problematic four wave mixing terms, denoted by reference numerals 1-10 (the amplitudes in all these plots are arbitrary and are intentionally displaced to make the terms easier to observe). Three of the four wave mixing terms fall directly on the shortest wavelength λ_(ref), two fall directly on the second channel λ₂, three fall directly on the third channel λ₃, and two fall directly on the fourth channel λ₄.

FIG. 8 shows an example in which the optical channels have irregular channel spacings (f₂−f_(ref)=250 GHz, f₃−f₂=200 GHz, f₄−f₃=150 GHz) and an optical channel filter bandwidth=50 GHz. The shortest channel wavelength is λ_(ref)=1310.000 nm. This plot shows the locations of 9 potentially problematic four wave mixing terms.

FIG. 9 shows another example in which the optical channels have irregular channel spacings (f₂−f_(ref)=250 GHz, f₃−f₂=200 GHz, f₄−f₃=350 GHz) and an optical channel filter bandwidth=50 GHz. The shortest channel wavelength is λ_(ref)=1310.000 nm. This plot shows the locations of 9 potentially problematic four wave mixing terms.

FIG. 10 shows yet another example in which the optical channels have irregular channel spacings (f₂−f_(ref)=150 GHz, f₃−f₂=100 GHz, f₄−f₃=200 GHz) and an optical channel filter bandwidth=50 GHz. Once again, the shortest channel wavelength is λ_(ref)=1310.000 nm. This plot shows the locations of 11 potentially problematic four wave mixing terms.

FIG. 11 show another example in which the optical channels have irregular channel spacings (f₂−f_(ref)=150 GHz, f₃−f₂=100 GHz, f₄−f₃=300 GHz) and an optical channel filter bandwidth=50 GHz. The shortest channel wavelength is λ_(ref)=1310.000 nm. This plot shows the locations of 10 potentially problematic four wave mixing terms.

With the various choices of irregular channel spacings shown in FIGS. 8-11 only the four wave mixing products generated by the unavoidable totally degenerate cases fall directly on the optical carrier frequencies. All of the other four wave mixing terms fall outside of the optical channel filter bandwidths in harmless locations posing no issues to the integrity of the transmitted information.

The previous discussion illustrated the techniques for reducing FWM using four channels. The following analysis will illustrate how these techniques can be generalized to any number of channels.

For purposes of this analysis a system with multiple optical channels can be expressed as:

$\begin{matrix} {{S(t)} = {\sum\limits_{i = 1}^{n}{\cos\left\lbrack {2{\pi\left( {f_{ref} + {\left( {i - 1} \right)\Delta\; f} + {\delta\; f_{i}}} \right)}t} \right\rbrack}}} & (1) \end{matrix}$ where i is the channel index number, n is the number of optical channels in the system, f_(ref) is the reference frequency (i.e. the lowest optical carrier frequency in the system), Δf is a standard frequency interval for the channel spacing, δf_(i) are the nonstandard frequency offsets, and t is time.

In DWDM (dense wavelength division multiplexed) optical communication systems Δf is typically either 50 GHz, 100 GHz, 200 GHz or even 400 GHz. In some cases channel separations as small as 25 GHz are also coming into use. The optical carrier frequency of the signal in the ith channel is given by: f _(i) =f _(ref)+(i−1)Δf+δf _(i)  (2) If we require that: δf₁=0  (3) then: f₁=f_(ref)  (4)

Furthermore if we also require the frequency spacing Δf and frequency offsets δf_(i) to be positive or zero: Δf≧0 and δf_(i)≧0  (5)

Then the optical carrier frequencies in (2) are ordered according to their index number: f₁=f_(ref)<f₂<f₃ . . . <f_(n)  (6)

The four wave mixing terms arise as a result of the third order nonlinearity in the fiber. Cubing (1) gives:

$\begin{matrix} {\left\lbrack {S(t)} \right\rbrack^{3} = {\sum\limits_{i = 1}^{n}{\sum\limits_{j = 1}^{n}{\sum\limits_{k = 1}^{n}{{\cos\left\lbrack {2\pi\; f_{i}t} \right\rbrack}{\cos\left\lbrack {2\pi\; f_{j}t} \right\rbrack}{\cos\left\lbrack {2\pi\; f_{k}t} \right\rbrack}}}}}} & (7) \end{matrix}$

Using appropriate trigonometric identities gives:

$\begin{matrix} {\left\lbrack {S(t)} \right\rbrack^{3} = {\frac{1}{4}{\sum\limits_{i = 1}^{n}{\sum\limits_{j = 1}^{n}{\sum\limits_{k = 1}^{n}\begin{Bmatrix} {{3{\cos\left\lbrack {2{\pi\left( {f_{i} + f_{j} - f_{k}} \right)}t} \right\rbrack}} +} \\ {\cos\left\lbrack {2{\pi\left( {f_{i} + f_{j} + f_{k}} \right)}t} \right\rbrack} \end{Bmatrix}}}}}} & (8) \end{matrix}$

The second term within the curly brackets in (8) represents the sum frequency terms which occur at frequencies in the vicinity of 3f_(ref) (approximately at ⅓ the reference wavelength). These terms are not problematic in a typical DWDM system. Therefore we may ignore this term and concentrate only on the first term within the curly bracket which produces spectral components that fall in the DWDM optical frequency band of interest. Using (2) this term becomes:

$\begin{matrix} {{S_{FWM}(t)} \equiv {\sum\limits_{i = 1}^{n}{\sum\limits_{j = 1}^{n}{\sum\limits_{k = 1}^{n}{\frac{3}{4}{\cos\left\lbrack {2{\pi\begin{pmatrix} {f_{ref} + {\left( {i + j - \left( {k + 1} \right)} \right)\Delta\; f} +} \\ {{\delta\; f_{i}} + {\delta\; f_{j}} - {\delta\; f_{k}}} \end{pmatrix}}t} \right\rbrack}}}}}} & (9) \end{matrix}$

The frequencies of the four wave mixing terms are given by: f _(i,j,k) ^(FWM) ≡f _(ref)+(i+j−(k+1))Δf+δf _(i) +δf _(j) −δf _(k)  (10a)

Notice that:

$\begin{matrix} \begin{matrix} {f_{i,i,i}^{FWM} = {f_{ref} + {\left( {i - 1} \right)\Delta\; f} + {\delta\; f_{i}}}} \\ {= f_{i}} \end{matrix} & \left( {10b} \right) \end{matrix}$

So that for the totally degenerate four wave mixing cases (i=j=k), the terms fall directly on the DWDM signals themselves. Those of skill in the art will appreciate that they are usually small enough in magnitude so as not to pose too much of a problem.

The indices i, j, and k range in value from 1 through n therefore: f_(1,1,n) ^(FWM)≦f_(i,j,k) ^(FWM)≦f_(n,n,1) ^(FWM)  (11a) Or: [f _(ref)+(1−n)Δf+2δf ₁ −δf _(n) ]≦f _(i,j,k) ^(FWM) ≦[f _(ref)+2(n−1)Δf+2δf _(n) −δf ₁]  (11b)

This may be further simplified using (3): [f _(ref)+(1−n)Δf−δf _(n) ]≦f _(i,j,k) ^(FWM) ≦[f _(ref)+2(n−1)Δf+2δf _(n)]  (11c)

The four wave mixing terms range over a spectrum that is three times the spectral occupancy of the DWDM signals themselves. f _(n,n,1) ^(FWM) −f _(1,1,n) ^(FWM)=3(f _(n) −f ₁)  (12)

Therefore, only a subset of these terms present a problem in that they may fall too close to and therefore interfere with any one of the original DWDM signals which are transmitting various types of information to the users of the system. We are thus interested in the four wave mixing terms which fall closest to any of the DWDM signal frequencies. We may examine the m^(th) signal. Using (2) and (10a) The frequency separation between the m^(th) DWDM signal and the four wave mixing terms are given by: ΔF _(i,j,k) ^(m) ≡f _(m) −f _(i,j,k) ^(FWM)=(m+k−i−j)Δf+δf _(m) +δf _(k) −δf _(i) −δf _(j)  (13) where 1≦m≦n. If i=j=k=m, clearly the frequency separation is zero which corresponds to the totally degenerate case also given by (10b). Other than the totally degenerate case the closest four wave mixing terms to f_(m) occur when the following condition is satisfied: m=i+j−k  (14)

Now notice that with the restriction (14), if i or j is equal to k then, it follows that j or i is equal to m. Hence we also have for the partially degenerate and totally degenerate cases: ΔF_(m,j,j) ^(m)=Δf_(i,m,i) ^(m)=ΔF_(m,m,m) ^(m)=0  (15)

Therefore, the totally degenerate and partially degenerate four wave mixing terms which are closest to the DWDM optical carrier frequencies fall directly on these frequencies and cannot easily be removed. Hence we concentrate on ensuring that the completely non-degenerate four wave mixing terms (i.e. i≠j≠k) fall outside of the optical filter passbands for each DWDM channel in the system.

Referring to (13) it reduces to: ΔF _(i,j,k) ^(m) =δf _(m) +δf _(k) −δf _(i) −δf _(j) with m=i+j−k, i≠j≠k  (16)

We can then eliminate the k index by using (14): ΔF _(i,j,(i+j−m)) ^(m) =δf _(m) +δf _((i+j−m)) −δf _(i) −δf _(j) with m≠i≠j  (17a)

Furthermore because the indices refer to the channel numbers according to (1) we must also impose the additional condition that: 1≦i+j−m≦n with m≠i≠j  (17b)

Typically in a channelized system, such as the WDM systems being discussed here, the optical signals are separated from one another at various switching sites within the communication network or at the optical receiver terminal sites so that the correct information can be transmitted to the proper consumer. This channel separation is usually accomplished by means of optical bandpass filters. Each filter has a given bandwidth of B and their respective bandwidths are centered on one of the n DWDM optical carrier frequencies f_(i). In order that no interference occur amongst the DWDM signals Each of the DWDM optical carrier frequencies must be separated from its nearest neighbors by at least B, the optical passband filter bandwidths. This can be accomplished by simply requiring that the standard frequency interval between the signals, Δf, be greater than B. Δf>B  (18)

The frequency offsets of the four wave mixing terms, ΔF_(i,j,k) ^(m), are determined with respect to the m^(th) optical carrier which, in-turn, is located at the center of its respective optical passband filter response. Therefore, in order to ensure that the generated four wave mixing terms do not fall within the passbands of the optical filters the magnitudes of the ΔF_(i,j,k) ^(m) must be greater than B/2.

$\begin{matrix} {{{\Delta\; F_{i,j,k}^{m}}} > \frac{B}{2}} & \left( {19a} \right) \end{matrix}$

Using (16) this becomes:

$\begin{matrix} {{{{\delta\; f_{m}} + {\delta\; f_{k}} - {\delta\; f_{i}} - {\delta\; f_{j}}}} > \frac{B}{2}} & \left( {19b} \right) \end{matrix}$

Or using (17):

$\begin{matrix} {{{{{\delta\; f_{m}} + {\delta\; f_{({i + j - m})}} - {\delta\; f_{i}} - {\delta\; f_{j}}}} > \frac{B}{2}}{{where}\text{:}}{1 \leq {i + j - m} \leq n}{and}{m \neq i \neq j}} & \left( {19c} \right) \end{matrix}$

Clearly if the signals are equally spaced so that δf_(m)=δf_(k)=δf_(i)=δf_(j) then the inequities (19) cannot be true. Hence one condition (although not sufficient) that needs to be satisfied in order that the four wave mixing terms fall outside of the optical filter passbands is that the channels must not be equally spaced. If the frequency offsets, δf_(i), are rewritten as:

$\begin{matrix} {{\delta\; f_{i}} \equiv {a_{i}\frac{B}{2}}} & (20) \end{matrix}$

Where the α_(i) are positive (or zero) scaling factors. Then condition (19c) can be expressed as: |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1 where: 1≦i+j−m≦n and m≠i≠j  (21) Also from (3): α₁=0  (22)

So long as condition (21) is satisfied for each of the frequency offsets as defined by (20) then the non-degenerate four wave mixing products will fall outside of the optical filter bandwidths. Also note that because (19c) and (21) are symmetric in the indices i and j that we need only consider unique pairs of (i,j) when performing the analysis. That is, in so far as the placement of non-degenerate four wave mixing products is concerned (i,j) gives the same result as (j,i).

For example, in a four wavelength system the unique (i,j) index pairings are given by the 10 pairs:

(1, 1) (1, 2) (1, 3) (1, 4) (2, 2) (2, 3) (2, 4) (3, 3) (3, 4) (4, 4)

Using equation 17(b), for each value of m the allowed values of i+j are given in table 4 along with the corresponding pairs of indices (i,j).

TABLE 4 Scaling Parameter Index Values (All Cases) All Allowed Corresponding M i + j Values (i, j) 1 2, 3, 4, 5 (1, 1), (1, 2), (1, 3), (2, 2), (1, 4), (2, 3) 2 3, 4, 5, 6 (1, 2), (1, 3), (2, 2), (1, 4), (2, 3), (2, 4), (3, 3) 3 4, 5, 6, 7 (1, 3), (2, 2), (1, 4), (2, 3), (2, 4), (3, 3), (3, 4) 4 5, 6, 7, 8 (1, 4), (2, 3), (2, 4), (3, 3), (3, 4), (4, 4)

TABLE 5 Scaling Parameter Index Values (Totally Non-Degenerate Cases m ≠ i ≠ j) Non Degenerate Corresponding M i + j Values (i, j) 1 5 (2, 3) 2 4, 5 (1, 3), (1, 4) 3 5, 6 (1, 4), (2, 4) 4 5 (2, 3)

Table 4 shows, for a particular value of m, all of the allowed i+j values and the corresponding values of the i and j indices. Table 5 isolates only the totally non-degenerate cases. For the four channel system there are six non-degenerate four wave mixing terms that are distinct and potentially problematic. Additionally, the partially degenerate and totally degenerate terms will generate terms falling directly on the four DWDM channel frequencies. Hence, in a four channel system there are 6 problematic, and distinct, non-degenerate frequency terms in addition to the indistinct partially degenerate and totally degenerate terms which fall directly on the 4 optical carrier frequencies. This gives a total of 10 distinct four wave mixing frequencies which may fall within the bandwidths of the optical channel filters:

TABLE 6 Number of Distinct Problematic Four Wave Mixing Frequencies Indistinct Total Number Degenerate Distinct of Distinct On-Carrier Non-Degenerate Four Wave FWM Four Wave Mixing Mixing m Frequencies Frequencies Frequencies 1 1 1 2 2 1 2 3 3 1 2 3 4 1 1 2

Table 6 summarizes the number of distinct problematic four wave mixing frequencies which occur near each carrier of a four channel DWDM system.

Table 5 can be used with (21) to define the frequency offset scaling parameter value inequalities:

|a₄ − a₂ − a₃ | > 1 (1, 2, 3) (23a) |2a₂ − a₃ | > 1 (2, 1, 3) (23b) |a₂ + a₃ − a₄ | > 1 (2, 1, 4) (23c) |a₃ + a₂ − a₄ | > 1 (3, 1, 4) (23d) |2a₃ − a₂ − a₄ | > 1 (3, 2, 4) (23e) |a₄ − a₂ − a₃ | > 1 (4, 2, 3) (23f)

In (23) the value in parentheses next to each inequality corresponds to the values of the indices (m, i, j) used to generate the inequality. Also (22) was used to simplify the expressions. Notice that (23a) is the same as (23c), (23d) and (23f). Then three remaining independent equations that remain are:

|a₄ − a₂ − a₃ | > 1 (1, 2, 3) (23a) |2a₂ − a₃ | > 1 (2, 1, 3) (23b) |2a₃ − a₂ − a₄ | > 1 (3, 2, 4) (23e)

Any values of α₂, α₃, and α₄ which satisfy these three inequalities (along with α₁=0) will generate non-degenerate four wave mixing terms that fall outside of the optical channel filter bandwidths and will not interfere with the transmission of the information being transmitted by the DWDM optical carriers.

One such solution is: α₁=0, α₂=2, α₃=2, α₄=0  (24)

Then if the optical channel filters bandwidths are each 50 GHz wide the corresponding frequency offsets are using (20): δf₁=0, δf₂=50, δf₃=50, δf₄=0GHz  (25)

In fact, for α₂=α₃=2 inequalities (23a,b,e) will each be satisfied for any value of α₄ which satisfies either of: 0≦α₄≦1 or α₄>5  (26a)

Implying that the frequency offset δf₄ should satisfy either of: 0≦δf₄<25GHz or δf₄>125GHz  (26b)

The solutions presented in equations 26a and 26b for a system with 4 channels corresponds to the solutions depicted in FIGS. 5 a and 5 b, which were discussed above.

As a second example consider the three channel case. Using equation 17(b), all of the unique (i,j) index pairings are given by the 6 pairs:

(1, 1) (1, 2) (1, 3) (2, 2) (2, 3) (3, 3)

TABLE 7 Scaling Parameter Index Values (All Cases) All Allowed Corresponding m i + j Values (i, j) 1 2, 3, 4 (1, 1), (1, 2), (1, 3), (2, 2), 2 3, 4, 5 (1, 2), (1, 3), (2, 2), (2, 3) 3 4, 5, 6 (1, 3), (2, 2), (2, 3), (3, 3)

TABLE 8 Scaling Parameter Index Values (Totally Non-Degenerate Cases) Non Degenerate Corresponding m i + j Values (i, j) 1 NONE NONE 2 4 (1, 3) 3 NONE NONE

In this case only the m=2 channel has a non-degenerate term that can be placed outside of the channel filter bandwidth by satisfying the single inequality:

$\begin{matrix} \begin{matrix} {{{{2a_{2}} - a_{3}}} > 1} & \left( {2,1,3} \right) \end{matrix} & (27) \end{matrix}$

So that either: 0≦α₃<2α₂−1 or α₃>2α₂+1  (28)

So that if α₂=1 then: 0≦α₃<1 or α₃>3  (29a)

Hence the frequency offsets are then (for 50 GHz channel filter bandwidths): δf₁=0, δf₂=25 and 0≦δf ₃≦25 or δf₃≧75 GHz  (29b)

As a final example, consider the five channel case. All of the unique (i,j) index pairings are given by the 15 pairs:

(1, 1) (1, 2) (1, 3) (1, 4) (1, 5) (2, 2) (2, 3) (2, 4) (2, 5) (3, 3) (3, 4) (3, 5) (4, 4) (4, 5) (5, 5)

TABLE 9 Scaling Parameter Index Values (All Cases) All Allowed Corresponding m i + j Values (i, j) 1 2, 3, 4, 5, 6 (1, 1), (1, 2), (1, 3), (2, 2), (1, 4), (2, 3), (1, 5), (2, 4), (3, 3) 2 3, 4, 5, 6, 7 (1, 2), (1, 3), (2, 2), (1, 4), (2, 3), (1, 5), (2, 4), (3, 3) 3 4, 5, 6, 7, 8 (1, 3), (2, 2), (1, 4), (2, 3), (1, 5), (2, 4), (3, 3), (2, 5), (3, 4), (3, 5), (4, 4) 4 5, 6, 7, 8, 9 (1, 4), (2, 3), (1, 5), (2, 4), (3, 3), (2, 5), (3, 4), (3, 5), (4, 4), (4, 5) 5 6, 7, 8, 9, 10 (1, 5), (2, 4), (3, 3), (2, 5), (3, 4), (3, 5), (4, 4), (4, 5), (5, 5)

TABLE 10 Scaling Parameter Index Values (Totally Non-Degenerate Cases) Non Degenerate Corresponding m i + j Values (i, j) 1 5, 6 (2, 3), (2, 4) 2 4, 5, 6 (1, 3), (1, 4), (1, 5) 3 5, 6, 7 (1, 4), (1, 5), (2, 4), (2, 5) 4 5, 6, 7, 8 (2, 3), (1, 5), (2, 5), (3, 5) 5 6, 7 (2, 4), (3, 4)

TABLE 11 Number of Distinct Problematic Four Wave Mixing Frequencies Indistinct Total Number Degenerate Distinct of Distinct On-Carrier Non-Degenerate Four Wave FWM Four Wave Mixing Mixing m Frequencies Frequencies Frequencies 1 1 2 3 2 1 3 4 3 1 4 5 4 1 4 5 5 1 2 3

Table 11 shows that there are 15 distinct non-degenerate four wave mixing frequencies in addition to the 5 indistinct degenerate terms which fall on the DWDM carrier frequencies for a total of 20 problematic four wave mixing frequencies that arise in a five channel system.

Next using the information in table 10 along with (21) and (22) form the inequality expressions for the 15 nondegenerate cases:

$\begin{matrix} {{{a_{4} - a_{2} - a_{3}}} > 1} & \left( {1,2,3} \right) & \left( {30a} \right) \\ {{{a_{5} - a_{2} - a_{4}}} > 1} & \left( {1,2,4} \right) & \left( {30b} \right) \\ {{{{2a_{2}} - a_{3}}} > 1} & \left( {2,1,3} \right) & \left( {30c} \right) \\ {{{a_{2} + a_{3} - a_{4}}} > 1} & \left( {2,1,4} \right) & \left( {30d} \right) \\ {{{a_{2} + a_{4} - a_{5}}} > 1} & \left( {2,1,5} \right) & \left( {30e} \right) \\ {{{a_{3} + a_{2} - a_{4}}} > 1} & \left( {3,1,4} \right) & \left( {30f} \right) \\ {{{{2a_{3}} - a_{5}}} > 1} & \left( {3,1,5} \right) & \left( {30g} \right) \\ {{{{2a_{3}} - a_{2} - a_{4}}} > 1} & \left( {3,2,4} \right) & \left( {30h} \right) \\ {{{a_{3} + a_{4} - a_{2} - a_{5}}} > 1} & \left( {3,2,5} \right) & \left( {30i} \right) \\ {{{a_{4} - a_{2} - a_{3}}} > 1} & \left( {4,2,3} \right) & \left( {30j} \right) \\ {{{a_{4} + a_{2} - a_{5}}} > 1} & \left( {4,1,5} \right) & \left( {30k} \right) \\ {{{a_{4} + a_{3} - a_{2} - a_{5}}} > 1} & \left( {4,2,5} \right) & \left( {30l} \right) \\ {{{{2a_{4}} - a_{3} - a_{5}}} > 1} & \left( {4,3,5} \right) & \left( {30m} \right) \\ {{{a_{5} - a_{2} - a_{4}}} > 1} & \left( {5,2,4} \right) & \left( {30n} \right) \\ {{{a_{5} + a_{2} - a_{3} - a_{4}}} > 1} & \left( {5,3,4} \right) & \left( {30o} \right) \end{matrix}$

Keeping only the unique expressions leaves the following:

$\begin{matrix} {{{a_{4} - a_{2} - a_{3}}} > 1} & \left( {1,2,3} \right) & \left( {30a} \right) \end{matrix}$ $\begin{matrix} {{{a_{5} - a_{2} - a_{4}}} > 1} & \left( {1,2,4} \right) & \left( {30b} \right) \end{matrix}$ $\begin{matrix} {{{{2a_{2}} - a_{3}}} > 1} & \left( {2,1,3} \right) & \left( {30c} \right) \end{matrix}$ $\begin{matrix} {{{{2a_{3}} - a_{5}}} > 1} & \left( {3,1,5} \right) & \left( {30g} \right) \end{matrix}$ $\begin{matrix} {{{{2a_{3}} - a_{2} - a_{4}}} > 1} & \left( {3,2,4} \right) & \left( {30h} \right) \end{matrix}$ $\begin{matrix} {{{a_{3} + a_{4} - a_{2} - a_{5}}} > 1} & \left( {3,2,5} \right) & \left( {30i} \right) \end{matrix}$ $\begin{matrix} {{{{2a_{4}} - a_{3} - a_{5}}} > 1} & \left( {4,3,5} \right) & \left( {30m} \right) \end{matrix}$

One possible solution can be found by letting α₃=1 and α₅=0 then letting α₂=2, α₄=5: α₁=0, α₂=2, α₃=1, α₄=5, α₅=0  (31a)

So that for a system employing optical channel filters with 50 GHz bandwidths the corresponding frequency offsets are: δf₁=0, δf₂=50, δf₃=25, δf₄=125, δf₅=0  (31b)

FIG. 12 is a flowchart showing one particular procedure for implementing the techniques described above for appropriately selecting the channel spacings in a system having n optical carrier channels so that FWM is reduced.

The method begins in step 405 by sequentially indexing the number of optical carriers (channels) in the system from f₁ through f_(n) in order from the lowest to the highest approximate frequencies (i.e., longest to shortest wavelength). The indices i, and j denote any of the optical channels f₁ through f_(n) and the index m denotes a given one of the optical channels.

In step 410 the lowest optical frequency is chosen to have a fixed value denoted the reference frequency f_(ref).

In step 415 an initial, uniform, standardized frequency spacing Δf (without offsets) is chosen between all of the optical carriers.

In step 420, all the distinct ordered pairs of the channel indices (i,j) are formed. For a system with n channels there are (n²+n)/2 such unique pairs. This includes the “diagonal” i=j pairs.

In step 425 only those ordered pairs of indices (i,j) are selected for which the restriction of equation (17b) is satisfied for all of the allowed values of the sums of the indices i and j for the m^(th) channel. That is, for each m, select those ordered pairs (i,j) such that 1+m≦i+j≦n+m, where 1≦i≦n and 1≦j≦n (as in the examples presented in Tables 4, 7, and 9).

In step 430, for each value of m, select only those ordered pairs of indices (i,j) determined in step 425 that satisfy the totally non-degenerate condition m≠i≠j i.e., none of the three indices are the same. (as in the examples presented in Tables 5, 8 and 10).

Having determined all of the pertinent non-degenerate index triples (m,i,j) for each channel m in the preceding steps, in step 435 use the inequality relationships presented in equation 21, |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1, to form all of the inequality relationships for the system. In addition, choose α₁=0 (as in equation 22) to simplify the set of inequality relationships.

In step 440, reduce the set of inequalities from step 435 by retaining only the unique inequality relationships, disregarding any duplicates or their equivalents (i.e. the negative expressions).

In step 445, determine the remaining frequency offset scaling factors α₂, α₃, . . . , α_(n), which satisfy all of the inequality relationships set forth in step 440. Although not necessary, in some cases it may be convenient to require the frequency offset scaling factors to be positive integers or zero.

In step 450 calculate the frequency offsets δf_(i) by multiplying each of the scaling factors determined in step 445 by the minimum acceptable frequency separation between the optical carrier frequencies and the nearest four-wave mixing term. This value might simply be B/2 where B is the bandwidth of the optical channel filters used in the system.

Finally, in step 455, the optical carrier frequencies are determined in accordance with the equation f_(i)=f_(ref)+(i−1)Δf+δf_(i). Channels located at these frequencies relegate the non-degenerate four-wave mixing products to locations that will not interfere with the system performance.

The processes described above, including but not limited to those shown in FIG. 12, may be implemented in a general, multi-purpose or single purpose processor. Such a processor will execute instructions, either at the assembly, compiled or machine-level, to perform that process. Those instructions can be written by one of ordinary skill in the art following the description of FIG. 12 and stored or transmitted on a computer readable medium. The instructions may also be created using source code or any other known computer-aided design tool. A computer readable medium may be any medium capable of carrying those instructions and include a CD-ROM, DVD, magnetic or other optical disc, tape, silicon memory (e.g., removable, non-removable, volatile or non-volatile), packetized or non-packetized wireline or wireless transmission signals.

A method and apparatus has been described for reducing four-wave mixing in an optical transmitter that can be used in a variety of WDM and DWDM optical networks. In particular, precisely selected wavelengths are utilized which force the FWM terms that are generated to lie far enough away from the desired signals so that they will not corrupt the information being transmitted by these signals. The technique described herein is readily applicable to optical transmitters that generate optical channels having wavelengths located at or near the zero dispersion wavelength of the transmission path through which they travel. 

1. A method of transmitting a WDM optical signal, comprising: modulating a plurality of optical channels that are each located at a different wavelength from one another with a respective one of a plurality of information-bearing electrical signals, wherein channel spacings between adjacent ones of the plurality of optical channels are selected such that the optical channels give rise to selected FWM terms that lie outside a channel bandwidth of the optical channels; multiplexing each of the modulated optical channels to form a WDM optical signal; and forwarding the WDM optical signal onto an optical transmission path, wherein the channel spacings between each channel i and its adjacent channel is equal to a constant frequency offset plus a frequency offset δf_(i), wherein the frequency offset δf_(i) is ${{\delta\; f_{i}} \equiv {a_{i}\frac{B}{2}}},$ where B/2 is a minimum acceptable frequency separation and α_(i) is a scaling factor that satisfies the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1 where 1≦i+j−m≦n_m≠i≠j, n is the number of optical channels, and i and j are indices denoting any of the optical channels between channel 1 and channel n, which, along with a given one of the optical channels m, produce the selected FWM terms.
 2. The method of claim 1 wherein the selected FWM terms include completely non-degenerate FWM products.
 3. The method of claim 1 wherein the optical channels are centered about a wavelength of 1550 nm.
 4. The method of claim 1 wherein α₁=0.
 5. The method of claim 4 wherein the scaling factor α_(i) satisfies only unique ones of the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1.
 6. The method of claim 5 wherein all values of the scaling factor α_(i) are equal to zero or a positive integer.
 7. The method of claim 1 wherein B is equal to the channel bandwidth of the optical channels.
 8. The method of claim 1 wherein the optical transmission path is located in a HFC network.
 9. The method of claim 1 wherein the optical transmission path is located in a CATV transmission network.
 10. The method of claim 1 wherein the optical transmission path is located in a PON.
 11. The method of claim 1 wherein the optical channels are located at wavelengths at or near a zero dispersion wavelength of the transmission path.
 12. The method of claim 11 wherein the optical channels are centered about a wavelength of 1310 nm.
 13. At least one computer-readable medium encoded with instructions which, when executed by a processor, performs a method including: selecting an initial set of channel spacings between at least three optical channels to be employed in a WDM system; and adjusting the initial set of channel spacings so that the optical channels give rise to completely non-degenerate FWM products that lie outside a channel bandwidth of each of the optical channels, wherein the three optical channel comprise n optical channels denoted channels f1, f2, f3 . . . fn, from lowest to highest in frequency, and i and j are indices denoting any of the optical channels between channel 1 and channel n, which, along with a given one of the optical channels m, produce the selected FWM terms and further comprising: determining selected non-degenerate index triples (m,i,j) for each optical channel m; and adjusting the channel spacings such that a frequency separation between the m^(th) signal and the FWM products arising from each of the selected triples (m,i,j) is greater than a minimum acceptable frequency separation, wherein f1 is located at a predetermined fixed frequency f_(ref) and frequencies f2, f3 . . . fn are located at frequencies determined in accordance with f_(i)=f_(ref)+(i−1)Δf+δf_(i) where 2<i≦n, Δf is a fixed frequency interval, and δf_(i) is a frequency offset, wherein the frequency offset δf_(i) is ${{\delta\; f_{i}} \equiv {a_{i}\frac{B}{2}}},$ number B/2 is a minimum acceptable frequency separation and α_(i) is a scaling factor that satisfies the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1 where 1≦i+j+m≦n and m≠i≠j.
 14. The computer-readable medium of claim 13 wherein the selected non-degenerate index triples (m,i,j) are determined by forming all (n²+n)/2 distinct ordered pairs of channel indices (i,j) and selecting only those channel indices (i,j) that satisfy the relation 1+m≦i+j≦n+m and m≠i≠j.
 15. The computer-readable medium of claim 13 wherein α₁=0.
 16. The computer-readable medium of claim 15 wherein the scaling factor α_(i) satisfies only unique ones of the inequality conditions a |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1.
 17. The computer-readable medium of claim 16 wherein all values of the scaling factor α_(i) are equal to zero or a positive integer.
 18. The computer-readable medium of claim 13 wherein B is equal to the channel bandwidth of the optical channels.
 19. The computer-readable medium of claim 13 wherein the optical channels are located in a bandwidth centered around 1310 nm.
 20. The computer-readable medium of claim 13 wherein the optical channels are located in a bandwidth centered around 1550 nm.
 21. A WDM optical transmitter, comprising: a plurality of optical sources for generating optical channels located at different wavelengths; a plurality of optical modulators each having a first input for receiving a respective one of a plurality of information-bearing electrical signals and an output coupled to an input of a respective one of the plurality of optical sources to thereby provide a plurality of modulated optical channels, wherein channel spacings between adjacent ones of the plurality of optical channels are selected such that the optical channels give rise to selected FWM terms that lie outside a channel bandwidth of the optical channels; and a multiplexer coupled to the plurality of optical sources to receive and combine the modulated optical channels to produce a multiplexed optical signal wherein the channel spacings between each channel i and its adjacent channel is equal to a constant frequency offset plus a frequency offset δf_(i), wherein the frequency offset δf_(i) is ${{\delta\; f_{i}} \equiv {a_{i}\frac{B}{2}}},$ where B/2 is a minimum acceptable frequency separation and α_(i) is a scaling factor that satisfies the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1 where 1≦i+j−m≦n_m≠i≠j, is the number of optical channels, and i and j are indices denoting any of the optical channels between channel 1 and channel n, which, along with a given one of the optical channels m, produce the selected FWM terms.
 22. The WDM optical transmitter of claim 21 wherein the selected FWM terms include completely non-degenerate FWM products.
 23. The WDM optical transmitter of claim 21 wherein α₁=0.
 24. The WDM optical transmitter of claim 23 wherein the scaling factor α_(i) satisfies only unique ones of the inequality conditions |α_(m)+α_((i+j−m))−α_(i)−α_(j)|>1.
 25. The WDM optical transmitter of claim 24 wherein all values of the scaling factor α_(i) are equal to zero or a positive integer.
 26. The WDM optical transmitter of claim 21 wherein B is equal to the channel bandwidth of the optical channels.
 27. The WDM optical transmitter of claim 21 wherein the optical channels are located at wavelengths at or near a zero dispersion wavelength of the transmission path.
 28. The WDM optical transmitter of claim 21 wherein the optical channels are located in a bandwidth centered around 1310 nm. 